OPEN 3 ACCESS Freely available online 



PLOS 



ONE 



Association of Gln27Glu and Arg16Gly Polymorphisms in 
Beta2-Adrenergic Receptor Gene with Obesity 
Susceptibility: A Meta-Analysis 



Hongxiu Zhang 1 *, Jie Wu 1 ' 2 , Lipeng Yu 3 



1 Department of Obstetrics and Gynecology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, China, 2 State Key Laboratory of Reproductive Medicine, 
Department of Obstetrics and Gynecology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, China, 3 Department of Orthopaedics, The First Affiliated 
Hospital of Nanjing Medical University, Nanjing, China 



Abstract 

Background: The beta2-adrenergic receptor (ADRB2) gene polymorphism has been implicated in susceptibility to obesity, 
but study results are still controversial. 

Objective: The present meta-analysis is performed to determine whether there are any associations between the Gln27Glu 
(rs1042714) or the Arg16Gly (rs1042713) polymorphisms in ADRB2 and obesity susceptibility. 

Methods: 'The PubMed (1950-2014), Embase (1974-2014), and China National Knowledge Infrastructure (CNKI, 1994-2014) 
databases were searched using the search terms ("Beta2-adrenergic receptor", "P2-adrenergic receptor" or "ADRB2"), 
"polymorphism," and "obesity". Fixed- or random-effects pooled measures were determined on the bias of heterogeneity 
tests across studies. Publication bias was examined by Egger's test and the modified Begg's test. 

Results: Eighteen published articles were selected for meta-analysis. Overall analyses showed that rs1042714 (Gln27Glu) was 
associated with significantly increased obesity risk in the heterozygote model (Gln/Glu vs. Gln/GIn: OR: 1.16, 95% CI: 1.04- 
1.30, / 2 = 49%, P = 0.009) and the dominant model (Gln/Glu + Glu/Glu vs. Gln/GIn: OR: 1.2, 95% CI: 1.00-1.44, / 2 = 55%, 
P = 0.04), whereas no significant association was found in the other models for rs1042714. Also, no significant association 
was found between the rs1042713 (Arg16Gly) gene polymorphism and the risk of obesity in all genetic models. In addition, 
neither rs1042713 (Arg16Gly) nor rs1042714 (Gln27Glu) showed any significant association with obesity susceptibility when 
the population were stratified based on gender. 

Conclusion: Our meta-analysis revealed that the rs1042714 (Gln27Glu) polymorphism is associated with obesity 
susceptibility. However, our results do not support an association between rs1042713 (Arg16Gly) polymorphisms and 
obesity in the populations investigated. This conclusion warrants confirmation by more case-control and cohort studies. 
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Introduction 

Compelling evidence demonstrates that both obesity and other 
related traits have a significant genetic component [1], and that 
these phenotypes result from an interaction between the genetic 
background and environmental factors [2]. The beta2-adrenergic 
receptor gene (ADRB2), as a lipolytic receptor in human fat cells, 
is associated with lipid mobilization. The most common single 
nucleotide polymorphism (SNP) occurs at codon 16 (Argl6Gly; 
rs 1042 7 13) and codon 27 (Gln27Glu; rs 1042 7 14). By altering the 
amino acid sequence in the extracellular N-terminus of the 
ADRB2, the rsl042713 and rsl042714 allele mutations are 
believed to alter ADRB2 function [3]. 



A number of polymorphisms have been well-studied in ADRB2 
and obesity. However, individual reports regarding ADRB2 
polymorphisms with obesity have produced inconsistent results. 
For example, a previous study including 4, 193 Japanese subjects 
indicated that rs 10427 13 was not a major contributing factor for 
obesity in Japanese men [4]. However, another study [5] showed 
that beta2-adrenoceptor polymorphisms may contribute to the 
development of obesity through gene-environmental interactions. 
Large et al. [6] found that in Swedish women obesity was 
associated with rsl042714, but not with rsl042713. The 
conflicting results of such studies may be a result of statistical 
underpower from sample sizes that were too small to detect any 
relationship between ARDB2 and risk of obesity. Therefore, we 
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performed a meta-analysis of all published case-control or cohort 
studies to clarify the association of ADRB2 polymorphism with 
obesity susceptibility. 

Materials and Methods 

Publication search and inclusion criteria and exclusion 
criteria 

The first report of significance of ADRB2 was published in 1954 
[7], therefore we selected the starting date of 1950 (and last search 
date of April 19, 2014) for the article search in Pubmed, Embase, 
and the China National Knowledge Infrastructure. The search 
terms used were: "Beta2-adrenergic receptor", "P2-adrenergic 
receptor", "ADRB2", "polymorphism," and "obesity". No 
language restrictions were imposed. For articles with overlapping 
data, we selected the publication with the most extensive data 
available. 

To be included in the meta-analysis, the identified articles had 
to meet all the following criteria: a) evaluation of Gln27Glu 
(rsl042714) or Argl6Gly (rs 10427 13) polymorphism and obesity, 
b) inclusion of quantitative information on the estimated risk of 
ADRB2 Gln27Glu (rsl042714) and/or Argl6Gly (rsl042713) 
polymorphism for obesity, and c) inclusion of complete informa- 
tion about all genotype frequencies, d) used a case control or 
cohort or cross sectional design, randomization or blinding is not 
necessary; The exclusion criteria were as follows: a) papers not 
related to ARDB2 polymorphism and obesity research, b) review 
articles, commentaries, or unpublished reports, c) papers without 
usable data, and d) duplicate publications. 

Data extraction and quality assessment 

We followed the Meta-analysis Of Observational Studies in 
Epidemiology (MOOSE) guidelines for reporting meta-analysis of 
observational studies [8] . The following items from each individual 
study were extracted: the name of the first author, year of 
publication, number of patients, gender, country, ethnicity, body 
mass index (BMI) cut point, sample size of obesity and control 
groups. The first stage was a review of titles and/ or abstracts for all 
identified citations, followed by a second review stage of full text 
publications. Two reviewers (HXZ and LPY) independendy 
assessed the eligibility of studies, and the third investigator (JW) 
arbitrated any disagreements by discussion and consensus. If allele 
frequencies were not provided, they were calculated from the 
corresponding genotype distributions. For information not avail- 
able in the published paper, relevant data was obtained by 
contacting the corresponding authors. Two reviewers (HXZ and 
LPY) also assessed independently rated the methodological quality 
of every included study by the "Newcastle-Ottawa Quality 
Assessment Scale" (NOS)[9]. This scale contains nine items (1 
point for each) in three parts: selection (four items), comparability 
(two items) and exposure (three items). 

Some authors provided data only on subjects of one gender, 
some authors gave information on subjects of each gender, while 
others failed to report gender at all (gender not identified). The 
latter studies were included only in the group of both genders 
combined. Then we calculated the Hardy- Weinberg equilibrium 
for every study, both in the main group and in the gender-based 
subgroups. Eventually, 18 publications were enrolled in the main 
analyses, including 1 5 case control studies, one random study, one 
cross sectional study and one cohort study. Using these genotype 
comparisons, we pooled together the populations of both genders 
from all studies, and performed gender-based subgroup analyses 
that included all suitable studies. Nine publications were included 
for the ADRB2 Argl6Gly (rsl042713) gender-based groups and 



17 publications were enrolled for the Gln27Glu (rsl042714) 
gender-based groups. 

Statistical analysis 

Summary statistics were estimated in Review Manager 5.1 
software (RevMan 5.1, Copenhagen: The Nordic Cochrane 
Center, The Cochrane Collaboration, 2011). The association 
between the rs 10427 14 gene polymorphism and obesity was 
compared by the odds ratio (OR) with its 95% confidence intervals 
(CIs). The statistical significance of the summary OR was 
determined with the Z-test. Five comparisons were performed 
between the two groups: frequency of allele (Gin vs. Glu), 
heterozygote comparison (Gin/ Glu vs. Gin/ Gin), homozygote 
comparison (Glu/ Glu vs. Gin/ Gin), dominant model (Gin/ Glu + 
Glu/ Glu vs. Gln/Gln) and recessive model (Glu/Glu vs. Gin/ Gin 
+ Gln/Glu) of ADRB2 Gln27Glu. For the rsl042713 (Argl6Gly) 
polymorphism, we used the same strategy, by replacing Glu27 
with Argl6. Five comparisons were performed between two 
groups: frequency of allele (Arg vs. Gly), heterozygote (Arg/ Gly vs. 
Arg/Arg), homozygote (Gly/ Gly vs. Arg/ Arg), dominant model 
(Arg/Gly + Gly/Gly vs. Arg/Arg) and recessive model (Gly/ Gly 
vs. Arg/Arg + Arg/Gly) of ADRB2 Argl6Gly. 

In consideration of the possibility of heterogeneity among the 
studies, a statistical test for heterogeneity was examined by the 
Chi-square-based Ojtest, and the significance was fixed at the level 
P<0.05. The inconsistency index / was also calculated to evaluate 
the variation caused by the heterogeneity. A high value of f 
indicated a higher probability of the existence of heterogeneity. 
A random-effects model (DerSimonian and Laird method) was 
used if substantial heterogeneity was detected (Ojstatistic: P<0.10; 
7 2 >50%). Otherwise, a fixed-effect model was applied in the 
absence of between-study heterogeneity (Q;Statistic: P>0.10; 
f< 50%). The significance of the pooled OR was determined 
by the Z-test, and the significance was set at P<0.05. Fisher's exact 
test was used to assess the Hardy-Weinberg equilibrium, and the 
significance was set at P<0.05. Potential publication bias was 
estimated using a funnel plot. Egger's linear regression test was 
used to evaluate the funnel plot asymmetry on the natural 
logarithmic scale of the OR (P<0.05 was statistically significant). 
We also further investigated the rs 10427 14 and rs 10427 13 gene 
polymorphism with obesity stratifying the population based on 
gender. 

Results 

Study characteristics 

According to the search strategy, 34 published articles were 
identified for potential inclusion with full text obtained for ADRB2 
polymorphism and obesity. Three articles were excluded because 
two papers were correspondences [10,11] and another because it 
was a review [12]. 12 studies were excluded because groups were 
not divided by BMI [6,13,14,15,16,17,18,19,20,21,22,23] and one 
paper was excluded due to the fact that data correlating ADRB2 
with obesity were not available [24] . Thus, 18 studies 
[4,6,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40] met our 
inclusion criteria. Of these, 17 papers involving 9,995 subjects 
genotyped at rsl042714 (all except Angeli et al. [39]) and 10 
studies [4,6,25,26,28,30,31,33,39,40] including 7,322 subjects 
genotyped at rsl042713 were included. A flow chart of study 
selection is shown in Fig. 1. The distribution of genotypes in all 
included studies was consistent with the Hardy-Weinberg equilib- 
rium. The key characteristics of these articles are summarized in 
Table 1. 
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PRISMA Flow Diagram 



Total number of citations retrieved 
from electronic searches ((n=34)) 



Studies excluded for 
correspondences (n=2) and review 
(n=l) 



Full-text manuscripts retrieved for 
detailed evaluation (n=3 1) 



ADRB2 gene data not shown on the 
basis of BMI values (n=12) 



Trials appear to meet Inclusion 
criteria (n=19) 



Articles excluded for AD RB2 data of 
obese people and controls not 
available ( n=l ) 



Trials included in the meta-analysis 
and systematic review (n=18) 



Figure 1. Flow diagram of articles selection process for ADRB2 gene polymorphism and obesity risk. 

doi:1 0.1 371 /journal.pone.01 00489.g001 



Methodological quality of 18 studies including our meta- 
analysis. Overall, the methodological quality of the 18 studies 
was modest. In general, the mixture studies including 15 case 
control studies, one random study, one cross sectional study and 
one cohort study. The studies failed to protect against selection 
bias: the definition of obesity used for the study is not uniform. 
None of the studies used secure methods for ascertainment of 
exposure. The majority of the studies provided evidence on the 
reliability of methods for outcome assessment; however, only 
several studies explicitly stated that outcome assessment was blind 
to exposure status. Finally, only two publications [4,39] including 
in our meta-analysis clearly declared that no conflict with interest, 
others did not mention that (TableSl). Methodological quality of 
18 articles enrolled in our study presented in Table 2. 

Meta-analysis 

In the meta-analysis, 17 studies involved the rs 10427 14 
(Gln27Glu) gene polymorphism, and 10 studies met our criteria 
for the rsl042713 (Argl6Gly) polymorphism. The frequency of 
occurrence of ADRB2 Gln27Glu/Argl6Gly allele in the popula- 
tion and their distribution among various populations are 
presented in table 3. A random-effects model was used if 



substantial heterogeneity was detected (Q_-statistic: P<0.10; 
/ >50%) or a fixed-effect model was applied in the absence of 
between-study heterogeneity (Q-statistic: P>0.10; 7 2 <50%). 
Meta-analysis of rs 1 042 7 1 4 (Gln2 7 Glu) and rs 1 042 7 1 3 (Arg 1 6Gly) 
polymorphism on risk of obesity are shown in Table 4 and Table 5, 
respectively. As shown in Table 4 and Fig. 2, in the analysis of the 
rsl042714 (Gln27Glu) gene polymorphism, the heterozygote 
model in our current study exhibited a significant difference 
(Gln/Glu vs. Gln/Gln: OR: 1.16, 95% CI: 1.04-1.30, 7 2 = 49%, 
P = 0.009), indicating that risk of developing obesity with Gln/Glu 
heterozygotes was 1.16 times higher than those with Gln/Gln 
homozygotes. Additionally, the dominant model also exhibited a 
significant difference (Gln/Glu + Glu/Glu vs. Gln/Gln: OR: 1.2, 
95% CI: 1.00-1.44, 7 2 = 55%, P = 0.04), suggesting that risk of 
developing obesity with the Gln/Glu plus Glu/Glu genotype was 
1 .2 times higher than with the Gin/ Gin homozygotes (Table 4 and 
Fig. 3). On the other hand, no significant correlation between 
obesity and Gln27Glu genetic variant in ADRB2 was found in the 
other three comparisons (Table 4). Our meta-analysis also showed 
that in all genetic models there was no significant association 
between Argl6Gly genetic variant in ADRB2 and the risk of 
obesity (Table 5). 
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Figure 2. Association between rs 1042714 (Gln27Glu) gene polymorphism and obesity risk under heterozygote model. (Gln/Glu vs. 
Gln/GIn: OR: 1.16, 95% CI: 1.04-1.30, / 2 = 49%, P = 0.009). 
doi:1 0.1 371 /journal.pone.01 00489.g002 



Evaluation of publication bias 

Publication bias was assessed by the funnel plot and Egger's test. 
The funnel plot (heterozygote Gln/Glu vs. Gln/GIn) showed no 
apparent evidence of publication bias (Fig. 4). There was also no 
significant difference in Egger's test for the allelic genetic model, 
which suggested that the probability of publication bias was low in 
the present meta-analysis (t = 0.84, P = 0.424). 



Subgroup analyses by gender in the current study 

Taking into account possible gender-specific roles in etiology 
[25], we conducted subgroup analyses by gender in the present 
study. After stratification for gender, no significant correlation of 
the rsl042714 (Gln27Glu) or rsl042713 (Argl6Gly) polymor- 
phisms to obesity was observed in any of the genetic models 
(Table 4 and Table 5, respectively). 
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Figure 3. Association between rsl 04271 4 (Gln27Glu) gene polymorphism and obesity risk under dominant model. (Gln/Glu + Glu/Glu 
vs. Gln/GIn: OR: 1.2, 95% CI: 1.00-1.44, I 2 = 55%, P = 0.04). 
doi:1 0.1 371 /journal.pone.01 00489.g003 
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Table 4. Meta-analysis of rs1042714 (Gln27Glu) polymorphism on risk of obesity. 
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OR 
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Abbreviations: OR, odds ratio; CI, confidence interval; / , Cochran's c-based Q-statistic test for assessing the heterogeneity (>50% indicates a substantial heterogeneity). 
doi:1 0.1 371 /journal.pone.01 00489.t004 



Discussion 

Several studies involved in ADRB2 polymorphisms and obesity 
have been published [4,5,6,41]. However, published results are 
controversial. For example, Pereira et al. reported that Argl6Gly 
(rsl042713) was not a major contributing factor for obesity in 
Japanese men; however, it is believedto have a significant 
association with obesity in Japanese women [4] . Large et al. [6] 
reported thatADRB2 gene polymorphisms were markedly associ- 
ated with obese Caucasian women. However, Echwald [34] and 
Oberkofler [31] found no association between ARDB2 and 
obesity. Specifically, a study performed by Oberkofler [31] 



concluded that the two polymorphisms of Gln27Glu (rsl042714) 
and the Argl6Gly (rsl042713) in the ARDB2 gene are not a major 
factor contributing to obesity in Austrian women. Moreover, 
Echwald et al. found the Glu27 polymorphism of ARDB2 gene is 
not associated with obesity in the population of Danish Caucasian 
men. However, Ehrenborg [42] found that individuals carrying 
the E27 allele and/or the G16 allele had significandy higher BMI, 
and that the E27 allele of the beta2-adrenoceptor gene is 
associated with slightly to moderately elevated BMI. 

In 2007, Gjesting and colleagues [23] conducted a case-control 
study and meta-analysis examining 7,808 white people for eirany 
association between ADRB2 polymorphisms and obesity. In their 



Table 5. Meta-analysis of rsl042713 (Arg16Gly) polymorphisms on risk of obesity. 



Comparisons rs1042713 (ArgieGly) polymorphism and obesity (lOstudies) OR 95%CI fC/o) P 



Allele frequency comparison (Arg vs. Gly) 


1.02 


0.95,1.10 


25 


0.52 


Gender-based subgroup analysis with men 


0.89 


0.74,1.08 


51 


0.24 


Gender-based subgroup analysis with women 


1.1 


0.99,1.23 


27 


0.08 


Heterozygote comparison (Arg/Gly vs. Arg/Arg) 


1.05 


0.93,1.19 


9 


0.39 


Gender-based subgroup analysis with men 


1.1 


0.91,1.31 


27 


0.32 


Gender-based subgroup analysis with women 


0.99 


0.82,1.20 


0 


0.92 


Homozygote comparison (Gly/Gly vs. Arg/Arg) 


0.95 


0.82,1.09 


24 


0.47 


Gender-based subgroup analysis with men 


1.1 


0.90,1.36 


44 


0.34 


Gender-based subgroup analysis with women 


0.71 


0.47,1.07 


53 


0.1 


Dominant model (Arg/Gly + Gly/Gly vs. Arg/Arg) 


1.02 


0.91,1.14 


15 


0.74 


Gender-based subgroup analysis with men 


1.1 


0.93,1.30 


46 


0.26 


Gender-based subgroup analysis with women 


0.92 


0.77,1.11 


0 


0.39 


Recessive model (Gly/Gly vs. Arg/Arg + Arg/Gly) 


0.92 


0.82,1.04 


8 


0.17 


Gender-based subgroup analysis with men 


1.04 


0.88,1.23 


0 


0.61 


Gender-based subgroup analysis with women 


0.67 


0.44,1.02 


71 


0.06 



Abbreviations: OR, odds ratio; CI, confidence interval; I 2 , Cochran's c-based Q-statistic test for assessing the heterogeneity (>50%indicates a substantial heterogeneity). 
doi:1 0.1 371 /journal.pone.01 00489.t005 
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Figure 4. Funnel plot for rs1042714 (Gln27Glu) gene polymorphism on heterozygote Gln/Glu vs. Gln/Gln. The funnel plot showed no 
apparent evidence of publication bias. 
doi:1 0.1 371 /journal.pone.01 00489.g004 



study, genotype distribution of ADRB2 Gln27Glu and Argl6Gly 
was provided according to diabetic people and non-diabetic 
subjects. However the data of ADRB2 data according to obesity 
and controls were not available. Hence the article was excluded 
from our present meta-analysis. Their analysis provided the data 
of BMI stratified according to ADRB2 Gln27Glu and Argl6Gly 
genotype and they did not find significant correlation between 
these beta2-adrenergic receptor variants (both)and obesity. Fur- 
thermore, they did not find that the quantitative trait analyses 
showed any effect of the variants on obesity-related traits. 

In 2008, Jalba and colleagues [43] performed a meta-analysis 
involving ADRB2 gene and obesity and conducted statistical 
analysis in three ways: the heterozygote comparison (Gin/ Glu vs. 
Gin/ Gin), the homozygote comparison (Glu/ Glu vs. Gin/ Gin) 
and the dominant model (Gin/ Glu + Glu/ Glu vs. Gin/ Gin). Their 
results suggested that rs 10427 14 might be a significant risk factor 
for obesity in Asians, Pacific Islanders, and American Indians, but 
not in Europeans. Also, the report showed that obesity might not 
be associated with rs 10427 13 at all. 

Since the association of ADRB2 polymorphisms with obesity is 
controversial, we conducted this meta-analysis based on all current 
available data on the relation between ADRB 2 polymorphism 
and obesity in 18 publications in order to clarify their relationship. 
Our findings suggest that there is a significant association between 
rsl042714 polymorphism of ADRB2 and obesity: OR = 0.86 for 
the allelic genetic model, OR =1.2 for the dominant genetic 
model, OR = 0.99 for the recessive genetic model, OR= 1.01 for 
the homozygote genetic model, OR= 1.16 for the heterozygote 
genetic model. 



In dominant comparison, the risk of developing obesity with 
Gln/Glu plus Glu/Glu genotype was 1.2 times higher than those 
with Gln/Gln homozygotes. However, our meta-analysis suggest- 
ed no significant correlation of the rs 10427 14 (Gln27Glu) 
polymorphism to obesity in the other three comparisons. Also, 
no significant correlation to obesity was found for the rs 10427 13 
polymorphism in all genetic models in the total population, as well 
as in gender-specific populations. This may be due to the 
divergence in genetic background. For example, the strength of 
the association between either rs 10427 13 or rsl042714 and 
obesity may be variable in different populations. Compared with 
the study of Jalba et al [43], we not only performed the following 
comparisons such as the heterozygote comparison (Gln/Glu vs. 
Gln/Gln), homozygote comparison (Glu/Glu vs. Gln/Gln) and 
dominant model (Gln/Glu + Glu/Glu vs. Gln/Gln), but we also 
conducted comparisons such as frequency of allele (Gin vs. Glu) 
and recessive model (Glu/Glu vs. Gln/Gln + Gln/Glu). In 
addition, we performed subgroup analyses by gender in the 
present study. Therefore, our results have stronger statistical 
power. Our meta-analysis suggest that Gln27Glu polymorphism of 
ADRB2 in heterozygote model showed a greater significance 
(p = 0.009) when compared to dominant model (p = 0.04). It may 
be a possible reflection of the linkage disequilibrium of genetic 
variability in codons 27 and gene-environment interaction in the 
etiology of obesity, since the mechanism of how Glu 27 can 
promote obesity is unknown at present. The results of our current 
meta-analysis support the conclusion that obesity susceptibility is 
associated with the Gln27Glu polymorphism of ADRB2 rather 
than the Argl6Gly polymorphism. This is different from previous 
results by Echwald [34] and Oberkofler [31] et al. Furthermore, 
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our findings are inconsistent with the results of Gjesting [23]. No 
association of ADRB2 (both) with obesity risk was observed in our 
meta-analysis upon gender stratification. However, our results 
support the conclusion fromjalba [43]. 

Though our study provides the most comprehensive and up-to- 
date meta-analysis regarding the association between ADRB2 
polymorphism with obesity, but also evaluates the methodological 
quality of 18 studies including our meta-analysis. Our work has 
some limitations. First, obesity is a complicated status involving 
complex interactions of genes, environment, and other factors, 
such as diet, lifestyle, diabetes mellitus, hypertension, total 
cholesterol, triglycerides, HDL and LDL. Though majority of 
publications included in our meta-analysis considered the factors 
described above such as diabetic factor, several articles did not 
consider the factors. See table 2. However, other susceptible 
factors were not able to be analyzed in the current study because 
insufficient data were provided from some of the original studies. 
Hence, misclassification bias is still possible. Moreover, only two 
studies [4,39] explicitiy stated that they have no conflict of interest, 
the others did not state whether they have conflict of interest or 
not. Therefore, possible conflicts of interests in studies enrolled in 
our meta-analyses may result in a bias. 

Second, different BMI values were used as a cut-off in enrolling 
studies. Asian populations require a lower BMI to indicate that an 
individual is at the same risk as a European, and this variation can 
be explained by WHO expert consultations [44]. Out of our 
included studies, BMI of less than 30 was used in four out of five 
Asian studies at codon 27 and all seven Asian studies at codon 16. 
The specific details of variations in BMI cut point are shown in 
Table 1 . It is not feasible to have the same BMI cut-off as obesity 
criteria across population from different geographical locations. 

Third, our results were based on an unadjusted estimate, a more 
precise analysis should be performed adjusted by age, smoking, 
and other factors. Lack of the original data of the enrolled 
publications limit our further evaluation of potential interactions 
such as gene-gene, gene-environmental factors, which may affect 
obesity risk. 
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